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Abstract 



> 
in 

\Q ■ It is shown that a weak modification of general relativity, in the lin- 

' earized approach, renders a spherically symmetric and stationary model 

, of the universe. This is due to the presence of a third mode of polarization 

■ in the linearized gravity in which a "curvature" energy term is present. 

' Such an energy can, in principle, be identified as the Dark Energy. The 

f^*) . model can also help to a better understanding of the framework of the 

Einstein- Vlasov system. 



PACS numbers 04.40.-b, 04.50.Kd, 04.40.Nr. 

The accelerated expansion of the Universe that is currently purported from 
observations of SNe la suggests that cosmological dynamics is dominated by a 
"new" substance of the universe constituents dubbed as Dark Energy, which is 
able to provide a large negative pressure to account for the late-time accelerate 
expansion. This is the standard picture, in which such a new ingredient is 
considered as a source of the right-hand- side of the field equations. It is posed 
that it should be some form of un-clustered non-zero vacuum energy which, 
together with the clustered Dark Matter, drives the global dynamics. This is the 
so-called "concordance model" (ACDM) which gives, in agreement with the data 
analysis of the observations of the Cosmic Microwave Background Radiation 
(CMBR), Lyman Limit Systems (LLS) and type la supernovae (SNe la), a 
good framework for understanding the currently observed Universe. However, 
the ACDM presents several shortcomings as the well known "coincidence" and 
"cosmological constant" problems [T] . 
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An alternative approach to explain the purported late-time acceleration of 
the universe is to change the left hand side of the field equations, and to inquire 
whether the observed cosmic dynamics can be achieved by extending general 
relativity [3M31EI- In this different context, it is not required to search candidates 
for Dark Energy and Dark Matter, which until to date, have not been found, 
but rather it claims that only the "observed" ingredients: curvature and baryon 
matter, have to be taken into account. Considering this point of view, one can 
posit that gravity is not scale- invariant [5]. In so doing, one allows for a room 
for alternative theories to be opened [HIEIIB]- In principle, interesting Dark 
Energy and Dark Matter models can be built by considering f(R) theories of 
gravity [5J[S] (here R is the Ricci curvature scalar). 

In this perspective, even the sensitive detectors of gravitational waves like 
bars and interferometers (i.e., those which are currently in operation and the 
ones which are in a phase of planning and proposal stages [TQl [IT], could, in 
principle, test the physical consistency of general relativity or of any other theory 
of gravitation. This is because in the context of Extended Theories of Gravity 
important differences with respect to general relativity show up after studying 
the linearized theory [JJ EH 031 US]- 

In this letter is shown that a weak modification of general relativity inspired 
by f(R) theories of gravity leads to a spherically symmetric and stationary model 
of the universe. (The analysis is performed in the context of the linearized the- 
ory) . Such a feature appears due to the presence of a third mode of polarization 
in the linearized gravity in which a "curvature" energy is present. This pecu- 
liar behavior may in some respect resemble the curvaton field dynamics. That 
is, it does not itself drive inflation; it merely generates curvature perturbations 
(which could be stationary) at late times after the inflaton field has decayed and 
the decay products have redshifted away, when the curvaton is the dominant 
component of the energy density. We recall that the curvaton field is a light 
scalar field during inflation whose quantum fluctuations produce the primordial 
density perturbations in a scenario for the origin of structure formation. It is 
posed that spatial variations in the curvaton density are then transferred to the 
cosmic background radiation density when the curvaton decays some time after 
inflation. See the following short list of references []Jj] . 

Such an energy can, in principle, be identified as the Dark Energy. The model 
can also help to have a better understanding of the Einstein- Vlasov system 



Equation ijTJ is a particular choice in f(R) theories of gravity [21 [3j [4j [5l [6], 



In cosmology the action ([T]) has been analysed in [21] in a rather different 
cosmological scenario as compared to the one analysed in this letter. In the 
limit e — * and /o = 1 it recovers the canonical form of the Einstein-Hilbert 
action of general relativity [21] [22], i-e. 



da nana ed!- 



Let us consider the action 




(1) 
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S = / (Txy/=gR + C m (2) 



Criticisms on f(R) theories of gravity arises from the fact that lots of such 
theories can be excluded by requirements of Cosmology and Solar System tests 
[2"5] . However, in the case of the action ([T]), the discrepancy with respect to 
the standard General Relativity is very weak, because e is a very small real 
parameter. Thus, the mentioned constraints could, in principle, be satisfied. In 
particular the authors of [23] found 

< e < 7.2* 1CT 19 . (3) 

We have also to emphasize that fundamental constrains can be renormalizcd in 
order to obtain fa = 1. 

Because we want to study interactions at cosmological scales, the linearized 
theory in vacuum, i.e. with C rn = 0, must be considered. Notice that it gives a 
better approximation than the Newtonian theory [28] and the importance of the 
linearized theory in a cosmological framework has also been recently emphasized 
by George Ellis [29). Therefore, we will analyse the pure curvature action 

v 4 ^. n7.t„ pi+e 



S = J d 4 x^g~f R L+£ . (4) 

Also notice that the theory arising from such an action has been recently 
linearized in [30] . but a review is needed for a better understanding of the 
theoretical framework. 

By varying the action Q with respect to the field equations arc obtained 
(through this paper the convention G = 1, c = 1 and h = 1 will be used) 



G» v = {l + ^ fQR A -lg^efoR 1+e + l(l+e)f R%^-g»M(l+e)foR £ }}- (5) 
By taking the trace of the field equations ([5]) one gets 

□ (l + £ )/ i? £ = ^y^/oi? 1+£ . (6) 

Then, by making the identifications [25l [30] 

$^(l+ £ )ioi? £ and % - ^/ i? 1+e (7) 

a Klein - Gordon equation for the effective $ scalar field is obtained. It can 
be written as 

— (8, 

To study gravitational waves, the linearized theory has to be analyzed, with 
a little perturbation of the background, which is assumed given by a near 
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Minkowskian background, i.e. a Minkowskian background plus <£> = <£>o (the 
Ricci scalar is assumed constant in the background) [13j [25] [30] . We also as- 
sume $o to be a minimum for the effective potential V: 

V ~ \a5<S> 2 ^^jr~ m 2 5<P, (9) 
2 d$ w 

and the constant m has mass dimension. 
Putting 

(10) 

$ = $0 + 5$. 

to first order in hnu and (5$, calling R^ vpa , i?^ and i? the linearized quan- 
tity which correspond to R pvpa , R pv and R, the linearized field equations are 
obtained [HHHISD]: 



Dh m = m 2 h m , 



where 



_ 5$ 

$0 

Then, from the second of eqs. pip , one can define the mass like [T31 1 



□ft m \35§ OSR 



5<f> V <5i? 



(11) 



(12) 



(13) 



Thus, as the mass is generated by variation of the Ricci scalar, we can say 
that, in a certain sense, it is generated by variation of spacetime curvature, 
re-obtaining the same result of [T3J [25] [30] . 

Note that, in the present case, the theory is suitable as the modification of 
General Relativity is very weak and in agreement with requirements of Cosmol- 
ogy and Solar System tests [33] . 

Rp,upa and eqs. (jlip are invariants under gauge transformations of the type 

nana [3D] 

hfj, v — > h^ u = — d( p e u ) 

(14) 

h E — > h! £ = h £ . 
Therefore, one can define [30] 



h pv = h pv - ^r\ pv + rj^he (15) 



and, considering the transform for the parameter as 
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□e„ = d" V. (16) 

On this basis one can choose a gauge similar to the Lorentz gauge that is used 
when studying electromagnetic waves. It reads: [3D] 

0" V = °- 

In this way, the field equations now are given by 

□V = ( 18 ) 
uh £ = (19) 

Solutions of eqs. JH]) and (19]) are plan waves [H ] [T3 ] 130] : 



hfj, v = A^lj?) exp(ip Q a; Q ) + c.c. (20) 
/i e = a(l?) exp(iq a x a ) + c.c. (21) 



where 



(22) 



Equation (|18p describes the dynamics of gravitational waves in standard 
general relativity [2TJ [22]. Equation ([20|) gives its solution. In parallel way, 
equations p^]) and ([2Tj) describe, respectively, the dynamics and the solution 
for the new mode (see also [T2l [T31 [50] ). 

It should be emphasized that the dispersion law for the modes of the "cur- 
vature" field h £ is not linear. Besides, the velocity of every "ordinary" wave 
mode h^v is the speed of light c, i.e., as it arises from general relativity. Rather, 
the second equation in (|22p) corresponds to the dispersion law for the mode 
h e . Because of this, the mass-energy field can be described like a wave-packet 
(T2l [T51 [3D]. Also, the group- velocity of a wave-packet of h e centered in in the 
momentum ~p* is 

«3 = -, (23) 

UJ 

which is exactly the velocity of a massive particle with mass-energy E and 
momentum j?. 

Therefore, from the second of eqs. ([2"2"]l and eq. ([2"5|) it is simple to express 
the group velocity as: 

v G = 1 ■ (24) 
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As one expects that the wave-packet possesses a constant speed, then it has 
to have an energy [12l EH [30] 



E = j0.-vl)u>. (25) 

On the other hand, the analysis can remain in the Lorenz gauge with trans- 
formations of the typede,, = 0; this gauge gives a condition of transversal effect 
for the ordinary part of the field: k^A^ = 0, but it does not guarantee the 
transversal effect of the total field h^ v . From eq. (p~5|) it becomes 

hfj, v = h^u - -Viiv + V^h e - (26) 

At this point, and considering that we are working in the massless case 
[121 H31 130] , this condition could be expressed as 

□e' 1 = 

(27) 

d^ = -± + h e , 

which provides the total transversal effect of the field. However, in the 
actual (massive) case this is impossible. In fact, by applying the D'alembertian 
operator to the second of eqs. ([27]) and using the field equations ([T8]) and ([T9l) 
one arrives to 

□e^ 1 = E 2 h £ , (28) 

which is in contrast with the first of eqs. (|27|) . In the same way it is possible 
to show that it does not exist any linear relation between the tensorial field 
and the "curvature" field h e . That is why a gauge in which /i M „ is purely spatial 
cannot be chosen (i.e. it cannot be given as h^o = 0, see eq. ((21))) ) . But the 
traceless condition to the field h^ v can be written as: [30) 

□e' 1 = 

(29) 

= -\. 

These equations imply 

d^h^ = 0. (30) 

In order to preserve the conditions d^h^ and h = one can use transfor- 
mations of the like 

□e^ 1 = 

(31) 

fye" = 0. 

By taking ~]f in the z direction, a gauge in which only An, A22, and A12 = 
A21 are different to zero can be chosen. The condition h = gives An = —^22- 
Now, substituting these equations in eq. l)2l))) . one obtains 
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/i„„(t, z) =A+(t- z)e i +) +A x (t~ z)e$ + h e (t - v G z)r,^. (32) 

The term A + (t — z)ejw + A x (t — z)e^t describes the two standard polar- 
izations of gravitational waves which arise from General Relativity, while the 
term h E (t — vqz)^^ is the massive field arising from the high order theory. In 
other words, the function R e of the Ricci scalar generates a third polarization 
state for gravitational waves which is not present in standard general relativity 
This third polarization has a "curvature" energy E. 

Now, a simple model of Universe will be proposed, in which the dynamic of 
the matter is described by the Vlasov equation [XT] [18l [19l [20] . However, the 
gravitational forces between the particles, viz., the galaxies, is now supposed to 
be mediated by the third mode of eq. ([3"2"]l after making the assumption that at 
cosmological scales such a mode becomes dominant (i.e. A + , A~ -C h e ) [2E]- In 
this way the "curvature" energy E can be identified as the Dark Energy of the 
Universe ~ 10~ 29 g/cm 3 [28l [29] . These two assumptions are exacltly the ones 
that concern the model of oscillating Universe in [55] . 

The model that we are going to discuss in this letter is parallel to the one 
introduced by Norstrom in [IS]. All the results will be obtained adapting the 
ideas introduced in [TTI ITS] \W \ 130] . 

In the hypothesis A + ,A~ <C h £ , the line element of our model will be the 
conformally flat one 

ds 2 = [1 + h £ (t, z)]{dt 2 - dz 2 - dx 2 - dy 2 ). (33) 

The possible using of a similar conformally flat line clement in a cosmological 
framework has been discussed in [28l [31] . 

To satisfy the condition demanding that the particles make up an ensemble 
with no collisions in the spacetime, the particle density must be a solution of 
the Vlasov equation 

*/ + £w-r^£<w = o. (34) 

Here V^ v represent the usual connections, / is the particle density and p° is 
determined by p a (a = 1, 2, 3) according to the relation 

WP" = - 1 (35) 

[T7l [T8l [191 [20] , which expresses the condition that the four momentum lies 
on the mass shell of the metric (greek indices run from to 3) [17] . 

We recall that, in general, the Vlasov- Poisson system is [17 ] [T8 ] IT9 ] [20] 

dtf + v- Vx/ - VxU ■ v J = 

AU = 4np (36) 

p(t,x) = J dvf(t,x,v), 
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where t denotes the time and x and v the position and the velocity of the 
galaxies. The function U = U(t, x) is the average Newtonian potential gener- 
ated by the galaxies. This system represents the non-relativistic kinetic model 
for an ensemble of particles with no collisions, which interacts through the grav- 
itational forces that they generate collectively [T71 [TBJ [TH1 ED] ■ Thus, one can use 
such a system to describe the motion of galaxies within the Universe, thought of 
as pointlike particles, when the relativistic effects are negligible [T71 fl8l fl9l [20] . 
In this approach, the function f(t,x,v) in the Vlasov- Poisson system is 
non-negative and gives the density on phase space of the galaxies within the 
Universe. 

The Vlasov equation (|34p implies that the function / is constant on the 
geodesic flow of the line element (|33p . The connection of the line element (|33[) 
are obtained from (note: as we are working in the linearized approach, in the 
following computations only terms up to first order in h e will be considered 
while high-order terms will be assumed equal to zero) 

= \{5^h e + 5«d v h e - —Lj-g^he). (37) 

In this way, the Vlasov equation in the spacetime defined by the line clement 
([33]) becomes 

*/ + £<W - + Tj^Jp]^/ = 0. (38) 

Now, let us recall that two quantities are important for the Vlasov equation 
in a curve spacetime |16| . The first is the current density 

N " = -J^V9P^f ( 39 ) 
and the second is the stress-energy tensor 

T^ = -J^^f,. (40) 

Here g is the usual determinant of the metric tensor that in the case of the 
line element (f33|) is given by 

.9=1 + 2h e . (41) 
Both of and arc divergence free (conservation of energy): 



(42) 

The mass shell conditions (|35|) can be rewritten as 

^(l + h £ )- 1 +5 ab p a p b . (43) 
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From the connections ([57]) , computing the Riemann tensor, Ricci tensor and 
Ricci scalar, the "effective" Einstein field equations 

Gfiv = TV, (44) 
can be obtained together with the "effective" Klein - Gordon equation 

□fts = -T, (45) 
where T = is the trace of the stress-energy tensor. 

To simplify the computations the analysis can be performed in a conformal 
frame. Thus, rescaling the stress-energy tensor in the form 

T?" = (1 + 3/i £ )T^, (46) 

one obtains 

T. = (1 + 2h e )T. (47) 

Note: in general, conformal transformations arc performed by rescaling the 
line-element like HZ] 

g a p = e*g a f3- (48) 
Here we choose the scalar field as being 

$ = h s , (49) 

which also implies 



in our linearized approach. 
Thus, equation (|45| becomes 



1 + h £ , (50) 



Uh e = — T*. (51) 

We note that the particle density is still defined on the mass shell of the 
starting line clement g a p- In order to remove even this last connection with the 
starting frame we rescale the momentum as 

Pt = (1 + yK (52) 
and define the particle density in the new conformal frame as 

f*(t,x,p«) = f(t,x, (1- y)p«). (53) 

Hence, we can write our adaptation of the Vlasov system in the following 
form 

nh e = (l+2h s ) J^-Mt,x,p*), (54) 
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P° = yl + 5abPipl (55) 

dtU + 4 9 -» /* - "4 [P*^^* + ^A]d p; /* = 0. (56) 
Pi P* 

Because we want to obtain spherical symmetry to remain in correspondence 
with astronomical observations, we can postulate the homogeneity and isotropy 
of the Universe. In this case the line-element ([55]) assumes the general form 

ds 2 = [l + h £ (t,r)]{dt 2 -dr 2 ). (57) 

where r is the radial coordinate. Thus, in spherical coordinates, equations 
(|5"^|) , (|5"5"|) and ([56]) can be written as 

-^F + re(|^) = (l + ^*.r), (58) 
/i(t,r)= / - 7 M=f(t,r, P ), (59) 

w+ 7!T7* / - [( ^7S7?s*' )pf 7n7'^* J ' w = °' <60) 

where the suffix * has been removed for the sake of simplicity, and we have 
denoted by p the vector p = (pi,P2,p3) with p 2 = |p| 2 , and also defined x for 
the vector = (xi, X2, X3). 

Then, wanting stationary states, following [28], we can call / the frequency 
of the "cosmological" gravitational wave ([52"]) . which has not to be confused 
with the particle density that we have previous denoted with /, and assume 
that / <§C Hq where Hq is the Hubble constant (i.e. the gravitational wave is 
"frozen" with respect the cosmological observations). In this case it is 4:h s = 0. 

Thus, the system of equations which defines the stationary solutions of 
cqs. ([58]) . ([59]) and ([60]) . for our model of Universe is 

— (^ = (1 + 2^), (61) 
r ar dr 

p(r)= / - 7 M=f(r,p), (62) 

pd r f--^-h e [{pr)p + r]d p f = 0, (63) 
r dr 



Conclusions 

It has been shown that a weak modification of general relativity, in the linearized 
approach produces a spherically symmetric and stationary model of the universe. 
This is due to the presence of a third mode of polarization in the linearized 
gravity in which a "curvature" energy is present. Such an energy can be, in 
principle, identified as the Dark Energy. The model can also help to a better 
understanding of the framework of the Einstein- Vlasov system. 
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